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Abstract. We determine the Hubble expansion and the general cosmic perturbation equa¬ 
tions for a general system consisting of self-conserved matter, pm, and self-conserved dark 
energy (DE), p£). While at the background level the two components are non-interacting, 
they do interact at the perturbations level. We show that the coupled system of matter and 
DE perturbations can be transformed into a single, third order, matter perturbation equa¬ 
tion, which reduces to the (derivative of the) standard one in the case that the DE is just 
a cosmological constant. As a nontrivial application we analyze a class of dynamical models 
whose DE density pd{H) consists of a constant term, Cq, and a series of powers of the Hubble 
rate. These models were previously analyzed from the point of view of dynamical vacuum 
models, but here we treat them as self-conserved DE models with a dynamical equation of 
state. We ht them to the wealth of expansion history and linear structure formation data 
and compare their ht quality with that of the concordance ACDM model. Those with (70 = 0 
include the so-called “entropic-force” and “QCD-ghost” DE models, as well as the pure linear 
model pd ~ 77, all of which appear strongly disfavored. The models with Co 7^ 0, in contrast, 
emerge as promising dynamical DE candidates whose phenomenological performance is highly 
competitive with the rigid A-term inherent to the ACDM. 
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1 Introduction 

It seems now beyond doubt that our Universe is in a state of accelerated expansion owing to 
some form of dark energy (DE) pervading all corners of the interstellar space. The evidence 
stems not only from the first and subsequent measurements of distant supernovae [1, 2] 
but also from the most recent analysis of the precision cosmological data by the Planck 
collaboration [3] . The physical cause for such positive acceleration is unknown but the simplest 
and most traditional possibility is to admit the presence of a tiny cosmological constant (CC) 
in Einstein’s field equations, A > 0. This canonical framework, the so-called “concordance” 
or ACDM model, seems to describe with good accuracy the available cosmological data [3] 
but, unfortunately, there is at present no special theoretical motivation for it. The CC 
term in Einstein’s equations is usually associated with the energy density carried by the 
vacuum through the parameter pA = A/(87rG) (in which G is the Newtonian coupling). The 
problem is that it is very difficult to reconcile its value (pA 10 GeV^) with the typical 
expectations in quantum field theory (QFT) and string theory. Such conflict between theory 
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and observation lies at the root of the famons and tongh “old CC problem” [4] and its variant 
the Cosmic Coincidence problem (see e.g. the reviews [5-8]). Both theoretical connndrnms 
stay at the forefront of fnndamental physics. Not snrprisingly there is plenty of motivation 
to search for alternative frameworks beyond the ACDM model. 

Needless to say the new proposed models, whatever it be their natnre, mnst be capable 
of endowing the A-term of a more satisfactory theoretical statns, and at the same time they 
shonld be able to maintain (or improve) the qnality of the hts to observational data. Different 
scenarios have been proposed in order to alleviate this sitnation, to wit: scalar helds, both 
qnintessence and phantom-like, modihed gravity theories, phenomenological decaying vacnnm 
models, holography scenarios, etc (cf. the previons review articles and references therein). In 
this work, we take the point of view that the dark energy density is a dynamical variable in 
Qnantnm Field Theory (QFT) in cnrved spacetime, as in snch context it shonld be possible to 
better tackle the basic CC problems In onr approach we do not acconnt for the dynamical 
character of the DE throngh specihc scalar held models or the like. We rather assnme the 
possibility that the effective form of the DE energy density in QFT in cnrved spacetime can be 
expressed as a generic power series of the Hnbble rate and its time derivatives. This possibility 
can be theoretically motivated from the renormalization gronp in cnrved spacetime, see e.g. 
[8, 22] and references therein. For the cnrrent Universe, however, the relevant terms can only 
be of order at most, whereas the important possible role played by the higher powers of 
H in the early Universe can be nsed to snccessfnlly implement inhation, see e.g. [9, 22-26]. 

The DE models nnder consideration in this work are closely related to those previonsly 
analyzed in detail in the comprehensive stndies of Ref. [15, 16], althongh with an important 
difference: now we consider that they describe a self-conserved dark energy, pd{H), with a 
dynamical eqnation of state (EoS) evolving itself with the expansion history: ojd = 

It means that pd{H) does not exchange energy with matter (which therefore remains also 
covariantly conserved). These assnmptions imply a completely different new class of DE 
scenarios which reqnires an independent cosmological analysis. We call them the “P-class” 
of dynamical DE models, to distingnish them from the vacnnm class (for which ud = —1 at 
any time of the cosmic expansion). We nndertake the task of examining in detail the P-class 
here. Of particnlar interest is the analysis of the effective EoS of the models in this class 
whose behavior near onr time conld explain the persistent phantom-like character of the DE 
withont invoking real phantom helds at all. 

In this paper we provide detailed considerations not only on the backgronnd cosmology 
of the new DE models bnt also on the pertnrbation eqnations and their implications on the 
strnctnre formation. The new pertnrbation eqnations are indeed formally different from the 
eqnations when the models are treated as dynamical vacnnm models in interaction with mat¬ 
ter. Of special signihcance is the formal proof that we provide according to which the conpled 

^For recent reviews of the idea of dynamical vacuum energy, see e.g. [8, 9]. See also previous works such 
as e.g. [10-17] and other recent and related studies such as e.g. [18-21]. 
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system of matter and DE perturbation equations for general models with self-conserved DE 
and matter components can be described by a single, third order, perturbation equation for 
the matter component. As a nontrivial application we subsequently solve (numerically) that 
equation for the P-class models and compare the results with the situation when the back¬ 
ground DE density is present but the corresponding DE perturbations are neglected. The 
main, and very practical, outcome of our work is that some of these dynamical DE models 
can provide a highly competitive fit to the overall cosmological data as compared to the per¬ 
formance of the concordance ACDM model - based on a rigid A-term. For the details of the 
specific cosmological observables used to fit our models to the data and the procedures of 
the statistical analysis we provide a short summary in Sect. 5.1. and refer the reader to the 
ample exposition previously put forward by some of us in [15, 16]. 

The layout of the paper is as follows. We present our dynamical DE models in Sect. 
2 and address their background solution and equation of state analysis in Section 3. The 
matter and dark energy perturbations are considered in Section 4. The confrontation with 
the background history and linear growth rate data is performed in Sect. 5. Finally, in Section 
6 we present our conclusions. 

2 Dynamical DE models: the P-class versus the vacuum class 

Let us consider a (spatially) flat FLRW universe. Einstein’s field equations read = 
8TTGTfj,i,, with G^i, the Einstein tensor and the total energy-momentum 

tensor involving matter and dark energy densities. We take both the matter and DE parts 
of the energy-momentum tensor in the form of a perfect fluid characterized by isotropic 
pressures and proper energy densities {pi,Pi)- One can then derive Friedmann’s equation 
and the pressure equation by taking the 00 and the ii components of Einstein’s equations, 
respectively: 

3H'^ = StT G {pd + PmPr) : (2-1) 

-|- 2H = —Stt G {pd + Pm, + Pr) j (^•2) 

where H = a/a is the Hubble function and the overdot denotes a derivative with respect to 
the cosmic time. As fluid components we consider cold matter, pm = 0, radiation, pr = Pr/3, 
and a general dark energy (DE) fluid with dynamical equation of state (EoS): pD = ojdPd 
{ uj£) ^ 0). If the matter and DE densities are separately conserved (in the local covariant 
form which we will indicate explicitly) we shall speak of self-conserved densities. Assuming 
also that there is no transfer of energy between non-relativistic matter and radiation (which 
is certainly the case in the epoch under study), the equation of local covariant conservation 
of the total energy density following from (2.1) and (2.2) can be split into three equations - 
reflecting the Bianchi identities of the Einstein tensor - as follows: 

Pm T 3HPm — 0 , 
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(2.3) 


Pr + AHpr = 0 , 


( 2 . 4 ) 


and 

Pd+ ^Hpui'^ + ojd) = ^ ■ (2.5) 

Similarly, we can e.g. check that from these local conservation laws and the pressure equation 
(2.2) we can reconstruct Friedmann’s equation (2.1). It follows that as independent set of 
equations we can take either the pair (2.1) and (2.2), or Friedman’s equation together with 
the local conservation laws. In the last case the pressure equation (2.2) is not independent. 
This will be our strategy in practice. 

Trading now the cosmic time derivative for the derivative with respect to the scale factor, 
through the simple relation d/dt = aH d/da, the above conservation laws can be easily solved: 

Pm{a) = p^^^a~^ , (2.6) 


Pr{a) = pi^'>a , (2.7) 

PD{a) = p^D exp i^-3 ^[1 + WD(a')]| > (2-8) 

where the densities pf'^ denote the respective current values. In the last case it is assumed 
that (jOd{o) has been computed. However, our method will be actually the opposite, we will 
compute the DE density Pd{o) first and then use (2.8) to compute the function ojd{o) - see 
Eq. (3.1) below. 

As it is transparent from (2.6) and (2.7), the pressureless (nonrelativistic) and the rela¬ 
tivistic matter energy densities are described by the standard ACDM laws, but the Universe’s 
evolution depends on the specific dynamical nature of the dark energy density po- The 
following basic DE models will be considered: 

VAI-. pd{H) = ^{Co + yH^) 

VA2 : pd{H) = + + ‘^aH^ (2.9) 

VA,: 

VCl : pd{H) = + uH^) 

VC2 : pd{H) = + \c^H ). (2.10) 

Notice that the constant, additive, parameter Cq has dimension 2 (i.e. mass squared) in 
natural units. We have introduced the dimensionful constant Hq (the value of the Hubble 
parameter at present) as a part of the linear term in H (for PCI) as in this way the free 
parameter e in front of it can be dimensionless. Similarly v and a are dimensionless parameters 
since they are the coefficients of and H, both of dimension 2. In the case of a we have 
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extracted an explicit factor of 2/3 for convenience. These free parameters will be fitted to 
the observational data. For all models we consider two free parameters at most. 

Models 2?A1 and PA3 in the list above are, of course, particular cases of 'DA2 cor¬ 
responding to a = 0 and = 0 respectively, but we have given them different labels for 
convenience and for further reference in the paper. Similarly, model PC2 is, too, a particular 
case of 2?A2 (for Cq = 0), but as we shall see it has some particular features that advice a 
separate study. Model PCI, in contrast, is not a particular case of any of the others. A 
particular realization of the PCI-model is the case u = 0 (i.e. the purely linear DE model in 
H), which will be denoted 

VH: pj,{H) = ^H. (2.11) 

This ostensibly simple model of the DE has been proposed in the literature on different 
accounts. It will be analyzed here in detail, along with the rest of the models, to ascertain 
its phenomenological viability (cf. Sect. 5). 

We remark that some of the above models are very similar to the ones we formerly 
called A1,A2 and C1,C2 in the comprehensive study of Ref. [15]. Formally the expressions 
for the DE densities are the same, the “only” difference being that in the previous reference 
they were treated as vacuum models (therefore with constant EoS, ujo = —1) in interaction 
with matter, whereas here the effective EoS is a function of the cosmological redshift, ud = 
u)d{z)^ with the additional feature that matter and DE are both conserved. The “P” in front 
of their names reminds us of the fact that these DE densities will be treated here in the 
fashion of dynamical DE models with a nontrivial EoS, the latter being determined by the 
equation of local covariant conservation of the DE density, namely Eq. (2.5). We will refer the 
cosmological DE models (2.9-2.11) solved under these specihc conditions as the “P-models”, 
or the models in the “P-class”, whereas we reserve the denomination of “dynamical vacuum 
models” when the same DE expressions are solved under the assumption that the EoS is 
uj£) = —1 at all times Let us also mention that the case of the pure linear DE model (2.11) 
was analyzed in detail in [16] from the point of view of a dynamical vacuum model, but here 
we will reassess its situation as a P-model^. 

The cosmological solution of the above P-models both at the background and perturba¬ 
tions levels turns out to be very different from their vacuum counterparts, as we shall show 
in this study. We refer the reader to Ref. [15] for the details on the vacuum part. 

^In this case, in order to fulfil the Bianchi identity, one has to assume that there is an interaction with 
matter [15] and/or that there is an additional dynamical component (as e.g. in the AXCDM model [27]), 
and/or that the gravitational coupling G is running [11, 21[. 

^In analogy with Ref. [15[ we could additionally have introduced the OBi models, namely the O-class 
analogous of the vacuum counterparts Bi defined there. The former are the model types obtained by replacing 
the H term of (2.9) with the linear term in H when Co yf 0. We will not address here the solution of the 
general models containing the linear term in H since it is not necessary at this point (cf. [15, 16] for details 
in the vacuum case). It will suffice for our purposes to study DCl and the pure linear model DH. 
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A few additional comments are in order before presenting the calculational details. On 
inspecting the various forms for the P-models indicated in (2.9-2.11), it may be questioned if 
all these possibilities are theoretically admissible. For example, the presence of a linear term 
oc Ff in some particular form of pd{H) - models DCl and, of course, DH- deserves some 
attention. Such term does not respect the general covariance of the effective action of QFT in 
curved spacetime [8] . The reason is that it involves only one time derivative with respect to 
the scale factor. In contrast, the terms and H involve two derivatives (d^ or d) and hence 
they can be consistent with covariance. From this point of view one expects that the terms 
oc and oc H are primary structures in a dynamical pD model, whereas oc Ff is not. Still, 
we cannot exclude a priori the presence of the linear term since it can be of phenomenological 
interest. For example, it could mimic bulk viscosity effects [18, 28, 29]. In addition, some 
of these models can be viewed also as holographic or entropic-force models [30], particularly 
the kind of models of PC2-type, which were proposed in Ref. [31]. Let us also mention the 
connection of the PH and PC-type models with the attempts to understand the DE from 
the point of view of QCD [32], and the so-called “QCD ghost dark energy” models and related 
ones [33-35]. In [12, 14] - see also [18] - these models were shown to be phenomenologically 
problematic. In fact, models PCI ,PC2 and PH present some phenomenological difficulties 
which we will elucidate here in the specific context of the P-class. In this respect let us 
emphasize that in order to identify the nature of these difficulties it is not enough to judge 
from the structure of the DE density, e.g. equations (2.9-2.11), as the potential problems 
may reside also in the assumed behavior of matter (e.g. whether matter is conserved or in 
interaction with the DE). That is why the list of pros and cons of the troublesome models 
examined in the previous references have to be carefully reassessed in the light of the new 
assumptions. We will accomplish this task here. As we will see, some of the old problems 
persist, while others become cured or softened, but new problems also appear. At the same 
time we will show that the only trouble-free models are indeed those in the large PA-subclass, 
i.e. the models possessing a well-defined ACDM limit. This was shown to be the case as 
dynamical vacuum models [15, 16] and it will be shown to be so, too, here as P-models. 


3 Cosmological background solutions 

The first task for us to undertake in order to analyze the above cosmological scenarios is to 
determine the background cosmological history. From the above equations it is possible, for all 
the models (2.9-2.11), to obtain a closed analytical form for the Hubble function, the energy 
densities in terms of the scale factor a or, equivalently, in terms of the redshift z = (1/a — 1), 
i.e. H{z). We use them to derive also the EoS and the deceleration parameters, which are 
very useful to investigate consistency with observational data. Thanks to the self-conservation 
of dark energy, the EoS can be extracted from the derivative of the DE density with respect 
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to the cosmic redshift 


(^oiz) = -1 + 


1 + z dpD{z) 


(3.1) 


2>pd{z) dz 

Similarly, the deceleration parameter emerges from the corresponding derivative of the Hubble 
rate, H: 

1 + z dH'^iz) 1 + z dE‘^{z) 

In the last expression we have used the normalized Hubble rate with respect to the current 
value, E{z) = H{z)/H q, a dimensionless quantity that will be useful throughout our analysis. 
The transition redshift point z = ztr from cosmic deceleration to acceleration can then be 
computed by solving the algebraic equation 


q{ztr)=^. (3.3) 

It is important to compute this transition point for the various models, as in some cases there 
may be signihcant deviations from the ACDM prediction. 

In the next section we systematically solve the background cosmologies for the PA- and 
PC-type models. The more complicated details concerning the solution at the perturbations 
level is presented right next. 


3.1 PA-models: background cosmology and equation of state analysis 

For the general PA2-model we have to solve the following differential equation: 

fju'i 

3(1 - u)H^ = + 3Co + aa — , (3.4) 

which follows from inserting the corresponding expression (2.9) of the DE density into Fried¬ 
mann’s equation and trading the cosmic time variable for the scale factor. We use also the fact 
that for the models under consideration the matter is locally self-conserved. The differential 
equation (3.4) is very different from the one obtained when the model is treated as a vacuum 
model in interaction with matter (cf. Ref. [15]) and therefore we expect that the PA-models 
should have a cosmic history different from the A-type ones. The constant Cq in (3.4) is hxed 
by imposing once more the current value of the DE density to be This yields 

- 1/ + ai^ + -I- 

where we have used the relations Hq = — (go + 1) Hq and 

go = -I- -I- (1 -I- , (3.6) 

with the EoS value of the DE today. Here the various stand as usual for the 
corresponding present day densities normalized with respect to the current critical density 


Co = Hi 
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= 3Hq/{8ttG). Upon integration of (3.4) we obtain the normalized Hnbble rate E{a) = 
H{a)/HQ, whose sqnare in this case reads: 


E^{a) = + 


Co 


Hiil- 






l — v + a 


(a ^-a^^) + 




(a (3.7) 


1 — 1 / + 4q;/3 


3/3 x 


where /3 = {1 — v)/a. We note the correct normalization iil(l) = 1. It is important to 
emphasize that we mnst have a > 0 for these models (in contrast to the sitnation with the 
A2-vacnnm type considered in [15]) since otherwise the term proportional to the derivative 
dH^/da on the r.h.s. of Eq. (3.4) conld become arbitrarily large and negative in the past, 
which wonld violate the non-negativity of in the corresponding l.h.s. of that expression. 
As a matter of fact, for the models with Co 7 ^ 0 (for which \u, a\ <C 1) we have /3 ;§> 1 since 
a cannot be negative. Thns, for the term appearing in (3.7) we obtain the nnll effective 
behavior ~ 0 for most of the cosmic history, namely nnless 2 ; is extremely 

close to zero. This observation does not apply for the models PC and PH, Eq. (2.9-2.11), 
since Co = 0 for them and hence the parameters a can no longer be simnltaneonsly small 
in absolnte valne. 

For PAl-type models we take the lateral limit a —>■ 0"*“ (i.e. we approach 0 from the 
right) in Eq. (3.7) (implying /3 —)> -|-oo), from which we can verify that all of the terms 
vanish since in this limit —)• 0 for a < 1. For a = 1 (the cnrrent time) these terms also 
cancel becanse the overall coefficient of is 0 in that limit. The Hnbble fnnction becomes, 
in this case, pretty mnch simpler: 


E\a) 


1 + 


1 - 


(a-3-l) + ^(a-^-l). 


(3.8) 


This resnlt for PAl can, of conrse, be obtained also from Eq. (3.4), which for a = 0 jnst 
becomes a simple algebraic eqnation. This shows the internal consistency of the obtained 
resnlts. For u = 0 we recover of conrse the ACDM resnlt. 

For illnstrative pnrposes let ns write down explicitly the evolntion of the DE density in 
the last case (i.e. for the PAl-model). Using the redshift variable we find: 


Pd{z) = 


rP 

Pd 


t^Pc 


1 - 


+ 




Pm (1 + + P^^'^ (1 + 2 :) 


(3.9) 


It is easily checked that it satisfies Pd{ 0) = p%, as it shonld, thanks to the cosmic snm rnle 
involving radiation: = 1. And of conrse it also boils down identically to 

for 1 / = 0 . 

Let ns come back to the more general model PA2. We neglect radiation at this point, 
as we want to focns now on featnres of the cnrrent Universe, snch as the eqnation of state of 
the DE near onr time. The normalized Hnbble fnnction (3.7) can then be cast in terms of 
the redshift as follows: 


EHz) = 


Co 


^ 1(1 - 


+ 


^ 7 ( 0 ) 


l — v + a 


{l + zf-7]{l + z) 
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Figure 1. DE densities normalized to their current value for the various dynamical DE models DA and 
DC under consideration. The behavior of DCl is different and has been plotted apart, together with the pure 
linear model DH. We note that in all plots (unless stated otherwise) we use the best fit values of Table 1 
corresponding to the barred quantities, i.e. those obtained when the structure formation data have also been 
taken into account in the fit. 


where 


T] = 


Co 


+ 


a 


( 0 ) 


Hq{1 — u) 1 — V + a 


- 1 . 


(3.11) 


The evolution of the DE density for the lDA2-models in the matter-dominated and current 
epoch can be computed with the help of the previous result, yielding 

P^cCo 


Pd{z) = 


Hiii- 


1 - zz -b a 


(l + z)3 -P°?7(1 + z)-3/3. 


(3.12) 


As can be checked, it satishes p_d(0) = p^{l — D^) = and it identically reduces to p^ for 
zz = a = 0, i.e. we recover in this limit the ACDM result. In the last part we use the fact that 
+ — 7 - 0 in the limit a —)• O'*" for any z > 0. One can verify that the PAl-type solutions 

(3.8) and (3.9) are particular cases of the last results in the matter-dominated and current 
epochs, as they should. Similarly, the DA3 -type solution is the particular case obtained 
from the above formulas for zz —)• 0. The numerical evolution of the DE energy densities for 
these models (normalized to the current value p^) is shown in Fig. 1 for the best ht values 
of Table 1 (see Sect. 5 for the details of the htting procedure leading to the results of that 
table). In Figures 2-3 we provide the corresponding behavior of the dark energy EoS and of 
the deceleration parameter for these models, which will be commented below in turn. 

From (3.12) and (3.1) we may derive the EoS for the general subclass of 'DA2 models. 
It can be expressed in the following compact form: 


u;d{z) 


P°cCo 
- v)pd{z) 




h(l + /3) 

Pd{z) 




(3.13) 


with Pd{z) given by (3.12). For Cq 7 ^ 0 the parameters zz and a are small (this is conhrmed 
by the htting values in Table 1) and therefore we can apply a similar argument to the limiting 
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Model 


n<"> 

UeE = l^-Ck 

i^eff 

o-g 

^8 

Xr/dof 

xVdof 

xVdof 

AIC 

AlC 

ACDM 

0.29+!!;“? 

0.286 ±0.007 


- 

0.815 

0.815 

569.21/592 

584.91/608 

584.38/608 

586.91 

586.38 

VAl 

0.286t“;“? 

0.281 ±0.005 

-0.024 ±0.018 

-0.028 ±0.016 

0.773 

0.770 

565.50/591 

573.02/607 

573.31/607 

577.02 

577.31 

VA2 

0.286=^0.011 

0.281 ±0.005 

-0.024 ±0.018 

-0.028 ±0.016 

0.772 

0.769 

565.57/591 

573.03/607 

573.40/607 

577.03 

577.40 

VA3 

0.287±0.011 

0.282 ±0.005 

U.UZO_g gjg 

-0.027 ±0.015 

0.777 

0.773 

565.63/591 

573.44/607 

573.47/607 

577.44 

577.47 

VCl 

0.286 ±0.014 

0.335 ±0.007 

-0.64 ±0.13 

—0.35 ± 0.05 

0.440 

0.735 

563.86/584 

880.74/600 

635.23/600 

884.74 

639.23 

VR 

0.242 ±0.008 

0.286 ± 0.005 


- 

0.513 

0.729 

639.85/585 

809.61/601 

677.11/601 

811.61 

679.11 

VC2 

0.285 ±0.013 

0.295 ±0.006 

J-VO_g 06 

1.02 ±0.01 

0.666 

0.752 

563.53/584 

594.13/600 

572.17/600 

598.13 

576.17 


Table 1. The best-fitting values for the various models and their statistical significance (x^-test and Akaike infor¬ 
mation criterion AIC[36], see Sect. 5). All quantities with a bar involve a fit to the total input data, i.e. the expansion 
history (BAO^i+BAOd^+SNIa) and CMB shift parameter data, as well as the linear growth data. Those without bar 
correspond to a fit in which we use all data but exclude the growth data points from the fitting procedure. The value Xr 
is the reduced which does not include the linear growth X/^g contribution. For models X^Al (resp. T>A3) i^eff = y 
(resp. —q); for T>K2 we have fixed a. = —v to break degeneracies (see text). For 7)0 and T>Yi models we have not used 
the BAOd^ and CMB data for the reasons explained in the text. In addition, for the T>Q models the given value of 
must be understood as the value of v since a. is not defined for 7)01 and becomes determined for 7)02 (see text). The 
quoted number of degrees of freedom (dof) is equal to the number of data points minus the number of independent 
fitting parameters. Details of the fitting observables are given in Sect. 5. 


situation just explained above to prove that the second term on the r.h.s. of (3.13) does 
not contribute significantly except in the origin. In practice, for any 2 ; > 0 (even for points 
very close to the origin) the effective EoS is actually given by the first term on the r.h.s. of 
(3.13). This is indeed the result that is connected by continuity with the EoS of the 'DAl- 
models in the limit a —)• O'*', as we shall show in a moment below. Therefore, for all practical 
situations related to redshift points around our current epoch we establish as effective EoS 
for the 2?A2-models the following expression: 


ujd{z) = 


1 + 




(l-u) 


Co 




ml-u+a 


(1 + 2 ) = 


(3.14) 


where by the same token we have disregarded the last term of Eq. (3.12). There is however 
one proviso related to the fact that Pd{z) in the denominator of the two terms in (3.13) could 
vanish, and in fact does vanish in our case for some particular redshift value (see below). This 
causes the presence of a vertical asymptote at a finite 2 point. In these cases, one would think 
of using Eq. (3.13) to better describe the behavior around the asymptote. In actual fact, not 
even this possibility affects in any significant way the practical use of Eq. (3.14), as we have 
checked. Being the parameters v and a small in absolute value we can expand the expression 
(3.14) linearly in them. The result can be cast in the following suggestive form for redshift 
points near our time (typically valid in the more accessible region 0 < 2 < 2): 


ujd{z) ~ -1 + 


Co 


z/eff (1 + zf 


-1 + 




1 - QO 


Z^eff (1 + zf 


(3.15) 


In this expression, the dimensionless quantity 


VeE = Tz -a 


(3.16) 
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is the basic fitting parameter. Using its value and that of from Table 1 we can evaluate 
the current EoS parameter of the VA2 model, Eq. (3.15): 

=a;D(0) ~-1.011, (3.17) 

which turns out to be remarkably close to the ACDM behavior. Let us also stress at this 
point that r'eg is indeed the effective fitting parameter for the entire set of DA-models at the 
background level in the matter-dominated and current epochs. It is a small parameter since 
\y,a\ <C 1 (owing to Cq 0). In particular, let us note that because is very close to 
— 1, the coefficient carrying Cq in the Hubble function (3.7) can be written in first order as 
in the limit —>• — 1. On the other hand the remaining u and a-dependence 

in (3 virtually disappears for the reasons discussed above. 

In the radiation epoch, however, the dependence on v and a is different than (3.16), as 
it is clear from the radiation term in (3.7). Also, in Sect. 5 we will see that this produces 
corrections to the transfer function which depend separately on u and a, and in these cases 
we must unavoidably fix some relation between these parameters. It has been indicated in 
the caption of Table 1, and more details are given in Sect. 5. 

There is one more, worth noticing, feature to stand out in connection to the numerical 
value (3.17) and the effective EoS function (3.15). They suggest that the dynamical DE 
models under study can provide, in principle, a reason for the quintessence and phantom¬ 
like character of the DE without necessarily using fundamental scalar fields. Indeed, we find 
that for v > a (i.e. t'eff > 0) the DA-models behave effectively as quintessence whereas for 
u < a (i.e. Ues < 0) they behave phantom-like near our time. In practice, after fitting the 
overall cosmological data, we have seen that the latter is the observed situation and also the 
predicted theoretical result. Recall that the current observational evidence on the dark energy 
EoS, for example from Planck results, leads to = —1.006 ± 0.045 [3]. This is perfectly 
compatible with the estimate (3.17) both quantitatively and qualitatively. This said, we are 
not suggesting that the current data on the EoS of the dark energy implies that the DE is 
phantom-like, as the accuracy around the central value is still insufficient. However, for years 
the central value of the measurement has shown some tilt into the region below —1, even 
during the long period of WMAP observations [37, 38]. Here we are merely saying that if 
such kind of measurement would be reinforced in the future, the general class of P-models 
encodes the ability for providing an explanation of the phantom character of the observed DE 
without need of invoking real phantom fields at all. Being r'eff = 0(10“^) rather small (and 
negative) the departure of —1 from below is very small, which is precisely the kind of result 
compatible with observations. To produce the plots of Fig. 2 (left) we have used both the 
exact expression (3.13) for the effective EoS of the PA2-model and the effective one (3.14), 
and have found no appreciable numerical differences. 

Let us now say some words on the presence of vertical asymptotes in Fig. 2. All models 
under study display these asymptotes, except PH (cf. next section). Such pole-like singularity 
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Cj(z) 




Figure 2. EoS function ujo{z) for the DE models VA (left) and for VCl and VC2 (right). The vertical 
asymptotes are located at the points where the DE density vanishes (compare with Fig. 1). The linear T>H 
model is seen not to present any asymptote (cf. right plot). 


is observed also in other contexts, for example in non-parametric reconstructions of the EoS 
function ojd{z) [39], in certain brane-model cosmologies [40] and in other situations, such as 
e.g. in mimicking quintessence and phantom DE through a variable A [13, 41]. In our case 
the pole-like feature is related to the fact that the denominator of (3.14) vanishes for some 
finite value of z, which is to be expected since t'eff < 0 (cf. Table 1). Physically this means 
that the DE density Pd{z) vanishes at these values of z (around 3.5 for models PA , and near 
1.5 or 2 for PCI depending on the fit options indicated in Table 1), as it is confirmed from 
the behavior of poi^) in Fig 1. As a consequence of this fact the EoS function (3.14) develops 
a singularity at this point. Notice, however, that the late-time expansion of the universe in 
the wide span 0 < z < 2 (possibly comprising all relevant supernovae data) is free from these 
exotic behaviors, if using the most optimized fit values that include the structure formation 
data. Obviously the latter are not associated with inherent pathologies of the model since the 
values of the fundamental physical quantities, such as the energy densities, stay finite (e.g. 
pniz) simply vanishes at these “singular” points). Interestingly, the very existence of these 
points might carry valuable information, for if the DE would be described by the P-models 
and we could eventually explore the EoS behavior at high redshifts {z > 2) we should be able 
to pin down these features, which would manifest through an apparent flip from phantom-like 
behavior into quintessence-like one when observing points before and after, respectively, the 
one where the DE density vanishes (cf. Fig. 1). Observing such phenomenon could be a 
spectacular signature of these models^. 

^It should be stressed that the present situation is different from the EoS studies previously entertained 
in Refs. [13, 41], in the sense that in the latter the effective EoS was only a representation (the so-called “DE 
picture”) of the original vacuum model (which in turn was called the “CC-picture” [41]), whereas here the EoS 
under study stands for the “physical” one associated with the original D-model. Therefore, if these models 
are to provide a correct description of the DE we should be able to observe the exotic EoS patterns shown in 
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Figure 3. Deceleration parameter, q{z), and transition point from deceleration to acceleration for the 
various DA and DC models. Notice that the transition point q{ztr) = 0 from deceleration to acceleration 
changes signihcantly for DC 1-models depending on whether we use barred or unbarred fitting parameters in 
Table 1. In contrast, DA and DC2-models have a more similar transition redshift ztr which is not far away 
from that of the ACDM. 


The EoS for the simpler subclass of PA 1-models can now be obtained from the limit 
a —)• 0^ of Eq. (3.13). In this limit the second term on the r.h.s exactly cancels for any 
z > 0 and we are left only with the first term. The result is simply Eq. (3.14) for a = 0. In 
particular, we can see that for z —>■ 0 we obtain a prediction for the current EoS value for this 
model: 


, ,(o) _ 

LOd — 


1 - - 1 / 


where the first expression is exact and the second is valid for |i^| <C 1. The latter is seen to be 
consistent with the limit z —)> 0 of Eq. (3.15). Since in this case the fitting values of u and i/gfr 
in Table 1 are the same, we retrieve the numerical result (3.17) also for PAl. Interestingly, 
we can explicitly verify that the result (3.18), which we have first obtained by taking the limit 
a —)• 0^ in Eq. (3.13), can also be worked out from (3.1) using the specific Eq. (3.9) of the 
PAl-model. 

The deceleration parameter for PA-type models reads: 


Zi —1 + u 




1 


(3.18) 


q{z) 


-1 + 


3Q 


( 0 ) 


2E^{z) — V + a 


(1 + z)=^ + 3/3r/(1 + z)-3^ 


(3.19) 


Solving Eq. (3.3) in this case we may find the transition redshift for a general PA2-type 
model. Once more we neglect the last term of Eq.(3.19) since it gives a negligible correction, 
and we arrive at 


Ztr = 


2(1 -b a - iy)[l - Pm"' - v + a{l + 


>(0) 

I 

Pm^(l - v) 


1 1/3 


- 1 


Fig. 2, much in the same way as in the alternative frameworks proposed in Refs. [39, 40]. 


(3.20) 
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The corresponding result for DAI is obtained by setting a = 0 in the above expression. Nu¬ 
merically, the deceleration and EoS parameters at the current time for DAI read respectively 
as follows: = —0.590 and ztr = 0.745. For DA3 models, the results are = —0.589 

and Ztr = 0.742. For v = a = 0 the formula (3.20) naturally returns the ACDM result: 


Ztr — 


f2^\ 



(3.21) 


whose numerical value for the htting parameters in Table 1 is ztr — 0.709. The plot of q{z) 
for the DA models is depicted in Fig. 3 (left), where we can also read the transition point 

Ztr- 


3.2 Models DCl, DH and DC2 


Let us start with DCl. The normalized Hubble rate E{z) = H[z)/Hq can easily be found 
starting from Friedman’s equation upon inserting in it the corresponding DE density from 
Eq. (2.10): 


E{z) 


e + T,{z) 
2 ( 1 -^) ’ 


(3.22) 


with 

S(z) = ^^62 + 4(1 - i/)[0^°^ (1 + z)4 + (1 + z)3] (3.23) 

The parameters in the above relation are constrained to satisfy 1 — = e -|- -j- so as 

to insure that E(0) = 1. Equivalently, e = — i/. Thus we can take e.g. i/ as the single free 

model-parameter, given the values of the ordinary parameters and For this model 

indeed, i/gfi in Table 1 is meant to be Furthermore, the same constraint shows that e and z/ 
cannot be both small parameters (in absolute value) in the case of the DCl models, i.e. they 
do not satisfy |e, z^| <C 1, in contrast to DA models, the reason being that the DCl models 
do not have a well-dehned ACDM limit for any value of e and z/. As a result one of the two 
parameter can be of order 0(1). For example, the ht in Table 1 for the barred quantities 
indicates that u ~ —0.35 and — 0.335, and therefore since the radiation component is 
negligible at present we obtain e ~ 1.02. For the case when the structure formation data 
are not used in the ht (unbarred parameters) we have ~ —0.64, which is much larger in 
absolute value, and then e is also larger: e ~ 1.35. 

The corresponding expression for the DE density of this model reads 


Pd(z) = pI [eE{z) + uE\z)\ , (3.24) 

where E(z) is given by (3.22). At 2 : = 0 the above function correctly renders the DE density 
now: Pd{ 0) = p^{€ + 1 ^) = p^ (1 — — D()) = Pc^% = P% upon using the mentioned 

constraint between the parameters of the model. For considerations in the matter-dominated 
and current epochs we can ignore of course the radiation component. A plot of the function 
(3.24), normalized to its current value is shown in Fig. 1 (right). In the same hgure we 
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plot the case of PCI when u = 0, i.e. the linear model PH, Eq. (2.11). This model has no 
free parameter, apart from since the above mentioned constraint enforces the relation 
e = 1 — matter-dominated epoch. Therefore, 

E{z) = - + S( 2 :)|^^q) = - + Y^py -I- 4P0j (1 -I- . (3.25) 

The corresponding DE density is: 

Pd{z) = Pc Pd E{z) = -p^D ^^y + Y^Py -I-dPj’j (1-|- z)^^ . (3.26) 


Notice the correct normalization /9d(0) = p^ since at z = 0 the argument in the square root 
becomes (2 — P^j)^ after using the sum rule P()j -|- P^^ = 1. It is interesting to remark at 
this point that the behavior of the PCI and PH models is rather different from their vacuum 
counterparts in interaction with matter. Let us e.g. focus on the linear vacuum model, i.e. 
the model (2.11) when the EoS is ojd = —1 and interacting with matter. As it is shown in 
Ref. [16], the corresponding DE and matter densities are 


Pa{z) —Pa |l + Pm (1 + ^ } 

Pm{z) = P^ p[(j -I- (1 — p[(j)(l -I-5;) (1-1- 2 )" 


(3.27) 

(3.28) 


In the matter density formula we note an extra factor of P[(j in front of the term (1 -|- z)^ 
for large enough z, as it was already pointed out in [16]. Thus, the deviation with respect 
to the ACDM behavior becomes increasingly large when we explore deeper our past. This 
is in contrast to its PH counterpart since matter is conserved for this model and hence it 
stays with the standard law Pm{z) = p^{^ + z)^. On the other hand, if we compare the 
dynamical DE densities (3.26) and (3.27), both deviate with respect to the ACDM linearly 
with z near our time (as can be seen by expanding these expressions for low z), but in in 
the PH case the relative deviation is a factor 1 — 2/(1 -|- P[(j) > 50% larger. This fact is 
actually determinant to explain the inability of the PH-model to correctly describe the low 
z data, and it reflects in a bad quality ht in Table 1. We conclude that the linear model 
(2.11) is essentially inefficient for a correct description of the cosmological data, both as a 
vacuum model and as a P-model. A similar situation occurs for the PCI and Cl models (cf. 
Table 1 and [16]). This is a clear sign that the P-models generally depart signihcantly from 
their vacuum analogs and both of them also with respect to the ACDM. In Sect. 5 we will 
retake this important issue in more detail, as these models have been repeatedly called for 
in the literature from different points of view and it may be appropriate to further discuss 
the reason for their delicate phenomenological status, which is made quite evident from the 
statistical analysis presented in the last four columns of Table 1. As we will see, it is not only 
caused by their background cosmological behavior but also by their troublesome prediction 
concerning structure formation. 
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Let us now address the EoS analysis of these models. Using (3.1) and (3.24) we are lead 
to the following expression for the dynamical EoS of the PCI and PH-models: 

Um^(l + z)^[e + 2i'E{z)] 


^d{z) — — 1 + 


E{z) [e + vE{z)] [2(1 - y)E{z) - e] ’ 


(3.29) 


As we are interested in the behavior of ujd{z) and q{z) in the low-redshift period (i.e. near 
our time), we have neglected the radiation contribution when we compute them. The corre¬ 
sponding EoS plot for PCI and PH are shown in Fig. 2 (right). In the case of PCI {v ^ 0), 
it displays an asymptote at large z near 2. The asymptote appears because the fitted value 
of V is negative in Table 1, so the last term of the expression (3.24) takes over at sufficiently 
large z and vanishes near z = 2 (cf. Fig. 1, right). There is no asymptote for PH because 
the function (3.26) is monotonically increasing with z. Moreover, one can show analytically 
that for large z the EoS for this model tends to —1/2, as can be confirmed graphically by 
looking at Fig.2. 

Of particular importance is the present-day value of the EoS parameter for PCI, i.e. 

Working out the result from the previous equations one finds: 




( 0 ) 

D 


e 

(1 — ^^m^)(e + 2Hm^) 


(3.30) 


Notice that this is an exact result and we recall that e and v are not small parameters for DCl 
models. Substituting the best-fit values of the parameters in the above formula, according to 
the results shown in Table 1, the current value of the EoS parameter achieved for PCI reads 
= —0.906^. This value can be spotted in Fig. 2 (right) and deviates significantly from 
the expected, narrow, range around —1 [3]. This goes, of course, also in detriment of model 
PCI. As for the pure linear model PH (corresponding to = 0), the constraint given above 
enforces e = (when we neglect the radiation component) and therefore Eq. (3.30) yields 

= (3.31) 

for the best-fit value of collected in Table 1. Such EoS result is out of the typical range 
of current observations (even more than for PCI) and puts the PH model once more against 
the wall. Let us also note that another particular case of PCI models, namely the pure 
quadratic /9 _d ~ model (which is obtained for e = 0), would have = 0 according to 
Eq. (3.30), and therefore it is completely excluded. The pure quadratic model was used in 
the past motivated by holographic ideas. But unfortunately the simplest ideas of this sort 
are actually ruled out [42] and some other generalizations (which include PCI-type models) 
too [12, 14, 18]. 

The deceleration parameter for the PCI-model can be computed from (3.2) and (3.22), 
and we find 

®Here and hereafter the numerical estimates use always the barred quantities in Table 1, which are the 
most optimized ones since they are obtained from fitting all the expansion and structure formation data. 
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(3.32) 


^ ^ [2 {I - u) - eHoH] ' 

The best-fit values indicated in Table 1 lead to the following estimates for the current 
acceleration parameters for the general PCI and PH, respectively: = —0.404, = 

—0.333. These models are signihcantly less accelerated now than the ACDM (for which 
q^^^ = —0.571). The behavior of q{z) is plotted in Fig. 3 (right). As it is seen in this 
hgure, for the PCI-model the Universe has a transit from a decelerated phase (q > 0) to an 
accelerated one ((7 < 0 ), which occurs precisely at the following transition redshift: 


^tr — 


2e2 




1/3 


- 1 , 


(3.33) 


Once more we borrow the htting results from Table 1 and for PCI models we obtain ztr = 
0.658, whilst for PH we have ztr = 0.528. These values are clearly smaller than for the ACDM 
(ztr — 0.691), meaning that the transition is accomplished much more recently than in the 
standard case. 

Let us now hnally deal with PC2-models. This model is sometimes also related with 
the entropic-force formulations, see [31]. However it is not enough to give the structure of the 
DE to know if the model is phenomenologically allowed or excluded. Treated as a vacuum 
C2-model {ojd = —1) it was shown in [14, 15] to be excluded, but as a P-model it has some 
more chances to survive. Let us describe here the main traits of its background evolution and 
leave the issues of structure formation for the next section. The Hubble function for PC2 
follows from (3.7) for Cq = 0. In the matter-dominated and current epochs reads 


E\a) = -(a-3 - . 

1 — u + a 


(3.34) 


As before f5 = {1 — v)/a, but we should recall that in this case a and u cannot be both 
small in absolute value. This is confirmed by explicitly displaying the constraint Cq = 0 that 
they satisfy in the matter-dominated epoch, namely 1 — — V + a = —in which 

~ —1. The evolution of the DE density for these models simply follows from (3.12) by 
setting Cq = 0, and we hnd: 


Pd{z) = p\ 


qI - V + a - ^ .,^- 3/3 , „ 0 r )0 


1 — V + a 


(l + z)-^^ + p[(L!^ 


v — a 
1 — p + a 


(i + T 


(3.35) 


As can be checked it satisfies /9 _d(0 ) = p^{l — D[(j) = if identically reduces to 

for u = a = 1, which is the ACDM result. Another way to see it is that, in this same 
limit, we have /3 —)• 0 and the normalized Hubble function (3.34) becomes exactly the ACDM 
one. Ultimately the reason for this is the following: for a DE density of the form PC2 in 
(2.10), the second Friedmann Eq. (2.2) can be satished identically if we set ojd = —1 and 
a = /3 = 1, and at the same time neglect the radiation component. Needless to say, the 
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self-conserved DE equation (2.5) is also automatically satisfied, because it is not independent 
of the two Friedmann’s equations. None of the other models could fulhl these conditions and 
consequently a nontrivial EoS different from ujd = —1 is required for them. It is thus not 
surprising that the htting procedure singles out a parameter region for 'DC2 very close to the 
ACDM (cf. Table 1). In actual fact the htting values of the parameters do carry some small 
deviations from unity, as this helps to slightly improve the quality of the ht as compared to 
the strict ACDM. This in turn induces a nontrivial evolution of the EoS, which eventually 
departs from ud = —1 when we move to higher redshifts, so in practice PC2 mimics the 
ACDM only near our time and quickly deviates from it for z > 1 (cf. Fig. 2, right). 

In contrast to its vacuum counterpart - viz. the C2 model studied in [14, 15] - the 
'DC2-model does have a transition point from deceleration to acceleration. A straightforward 
calculation yields 


Ztr = 


(0) 


— UJ 


( 0 ) 


D 




(3 — 3z^ -|- 2a) 


3{l+/3) 


- 1 = 


{1 — v + a — D((j)(3 — 3zz -|- 2 q;) 


3(l-i/+a) 


- 1 


(3.36) 

As expected, for v = a = 1 (which implies w® = — 1 from the aforementioned constraint) 
we recover exactly the ACDM result (3.21). As for the effective EoS of this model, ujj:){z), it 
can be directly computed with the help of (3.1) and (3.35). In particular one can check that 
for z = 0 we obtain a;£)(0) = —(1 — — y + a)/(aD2)) = which is consistent with the 

constraint obeyed by the parameters of this model. For the numerical analysis presented in 
Table 1 we have hxed = — 1 for this model and we have htted and zz, and in this 
way a becomes determined: a = {u + — 1)/D[(j. For this reason we can use u as the only 

vacuum free parameter for 2?C2 , and indeed this is the value that is meant for y^s in Table 1. 
The plots for the densities, the EoS behavior, as well as for the deceleration parameter and 
the transition point from deceleration to acceleration for PC2 are included in Figures 1-3. 
The model works well at low 2 ; but we should warn the reader that it may present serious 
difficulties in describing the radiation epoch, and in fact this is the reason why we have used 
only the low 2 ; observables in its ht. This fact puts us on guard concerning the eventual 
viability of 2?C2 . We will come back to this important point in Sect. 5.5. 


4 Linear structure formation with self-conserved DE and matter 

In this section we deal with the cosmic perturbations for linear structure formation, and we 
commence with a study of the set of perturbation equations for a general system of self- 
conserved dynamical dark energy and matter. Subsequently we specialize our results for the 
concrete P-models whose background behavior has been elucidated in the previous sections, 
all of them falling in that category. 
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4.1 Linear perturbations for matter and dark energy 


Let us consider a general system in which the various components (matter and DE) are 
covariantly self-conserved and the gravitational coupling G remains constant throughout the 
cosmic history. In the study of linear structure for this system we take into account both 
the matter perturbations, dpm, and the perturbations in the DE component, p^. Although 
the DE is not coupled to the matter sector at the background level owing to the assumption 
of separate covariant conservation, the two sectors develop some kind of interaction in the 
linear perturbation regime, which we are going to compute. We start by perturbing the 
(0,0) component of Einstein’s equations and the (0, i) components of the covariant energy 
conservation equation, i.e. = 0, using the synchronous gauge, i.e. ds^ = dt^ -|- 

{—a^Sij + hij)dx'^, and taking into account that all the components of the cosmic fluid are 
covariantly self-conserved. We obtain the following coupled system of diferential equations®: 

h 2Hh = 8 ttG [6pm “k dpf){l Slod) SpdSujd] (d-1) 


Pm 




+ 3H5pm + dpm — 0 


(4.2) 


Pd( 1 -k cod ) 



^ + ^H[{1 -k ujrdjdpD -k pd^uju] + dpD — 0 


(4.3) 


Pmidm 5140m) “k 0mPm — 0 


(4.4) 


|pd(1 + ^d){0d + 5H6d) + 6d [pd(1 + wd) -k pd^d]'^ 


a 

¥ 


{udSpd + PdSujd) = 0, (4.5) 


where On ^ divergence of the perturbed velocity for each component (matter 

and DE), h = -^ ^ stands for the wave number (recall that these equations are 

written in Fourier space [43]). 


4.2 Reduction to a single third-order differential equation for matter perturba¬ 
tions 

Notice that we have six unknowns and hve equations, so at this stage is not possible to solve 
the system in order to hnd each of the perturbed functions 6m, On, h, 5 ojd, SpD and 6pm- 
However, we will now show that under reasonable assumptions and working at sub-horizon 
scales we can eliminate the perturbations in the DE, i.e. 6pD, in favor of a single higher order 
equation for the matter part, 6pm- This is characteristic of the coupled systems of matter 
and DE perturbations for cosmologies with self-conserved dynamical DE and matter. For the 

®See e.g. Refs. [43-45] for a formulation of the perturbation equations in a system of several components, 
including the DE. Here we have introduced explicitly the perturbations in the dynamical EoS variable of the 
DE. 
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particular case A =const. the obtained third-order equation boils down to the (derivative of 
the) second order one of the ACDM, as we will show. Let us proceed step by step. Firstly, 
we make use of (2.3) and (4.4) so as to obtain 0m(a): 

0m{a) + 2H{a)9m{a) = 0 9m{a) = . (4.6) 

It follows that the matter velocity gradient decreases fast with the expansion. We shall adopt 
the conventional initial condition that is negligible or zero at present and hence 0^ — 0 
throughout the evolution. Since the scales relevant to the matter power spectrum remain 
always well below the horizon, i.e. k ^ H, we expect small DE perturbations at any time, 
i.e. the DE should naturally be smoother than matter, wherefore the velocity gradient of the 
DE perturbations is also naturally set to 0£) = 0. In this way we can conhne our study to 
the perturbation 6d in the DE sector and 5m in the matter sector. This simplihcation looks 
reasonable and it will allow us to solve the coupled system of matter and DE perturbations 
without introducing further assumptions and/or additional parameters. Under this setup it 
is clear that the terms that are not proportional to in (4.5) can be neglected and we hnd 
ojd^Pd + Pd^ujd = 0) which is tantamount to say that we are keeping the perturbations in 
the DE density but neglecting the perturbations in its pressure, similar as for matter. From 
(4.2) we can isolate h 

h= — (^H5pm + 5pm) (4.7) 

pm ^ ^ 

and differentiating this expression we arrive at: 

h = ^ (3H5pm + ^H5pm + 5pm^ - (sHSpm + 5pm^ ■ (4.8) 

Now we insert the last two equations in (4.1) and hnd: 

B{5pm, 5pm, Spm, H, H) = SttG [5pm + 5pd{1 + 3ujd) + 3pd5u}d] , (4.9) 

where we have dehned 

r- x" rrij^-^H5p m + 6H5p m + 25 p m Pm{QH5p m + 25 p m ) AH{m5p m + 5p m) 

B[5pm, 5pm, 5pm,B,H) = -—-— _____— -—+ ---— 

Pm Pm Pm 

(4.10) 

On the other hand from (4.7) and (4.3) we gather 

— —— (1 -|- LO£)){?)H5pm + 5pm) + 3iL[(1 -|- 0Jo)5pD + Pd^^d\ + ^Pd = 0 . (4-11) 

Pm 

We have three independent equations for the three unknowns: 5ojd-, ^Pd and 5pm- We hrst 
take 5ojd = —^d^Pd/Pd from the above mentioned equation and substitute it in (4.9) and 
(4.11). We hnd: 

B{5pmi ^Pm^ 5pmi H, FT) = 8ttG {5pf) 5pm) (4-12) 

and 

— ——(1 -|- ijj£)){2>H5pm + 5pm) + 3iF 5p£) -\- 5pD = 0. (4-13) 

Pm 
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At this point we can use the last two equations to get rid of the explicit 5pD dependence on 
the DE perturbations, and we obtain: 

“ + ^D){^H6pm + hpm) + ( — Spm\ + „ — Spm = 0 . (4-14) 

Pyyi yOTrCr J OTTCz 

This is the preliminary expression of the hnal differential equation for the linear density 
perturbations of the matter held, but in order to make it operative we must still do some 
algebra. Let us hrst rewrite the expression B, Eq. (4.10), as follows: 

B=^-PlIL{QH + 30H^) + —6pm + 2^-^, (4.15) 

pm Pm pm 

and compute its time derivative, 


B = ^{6H + 78HH + 90H^) + ^{78H^ + 22H) + 22H^ + —6pm , (4.16) 

Pm Pm pm pm 

where in the previous formulas we have made repeated use of the matter conservation (2.3). 
Introducing the last two expressions in (4.14), we obtain the hrst explicit form of the sought- 
for third order equation for 5pm- 


5pm + 96HH + 180H^ 
Pm SttG 


‘8H [pra + Pd(1 + ^d)]) 


25p m ^ 28H5p m 

SttGp m 8 'kGp m 


+ 


SPm 

Pm 


I2m^ + 22H 


-[{P m + Pni^ + wd)] 


= 0 . 


(4.17) 


A hnal simplihcation of Eq. (4.17) can still be performed. From the pair of Friedmann’s 
equations (2.1) and (2.2) we hud the relation 


H — —AirGlpm + P_d(1 + ^d)] ■ 


(4.18) 


Using it in (4.17) we hnally arrive at the desired form of the third order differential equation 
for the perturbations of the matter held: 

6pm + UH6pm + 3hp™(4i7 + 21H^) + 35pm{H + 30H^ + 17HH) = 0 . (4.19) 


It is convenient to reexpress it in terms the density contrast 5m = 5pm/Pm- After some 
algebra we hud the following rather compact form: 

(5 m + 7)H5m + 35m{H + 2H^) = 0. (4.20) 


This is actually the hnal form in terms of the cosmic time. However, to make contact with the 
observations it is highly advisable to rewrite the previous equation in terms of the scale factor, 
as it has a simple relation with the cosmic redshift. To this end we have to trade the time 
derivatives for the scale factor derivatives (indicating the latter by a prime), starting from 
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6m = ciHS'm and its subsequent second and third order derivatives. After a simple calculation 
we obtain: 


C + C[- + 


3H' 

~ir 


, 12 12H' F'2 jj/n 

+ I ^ H- rw + UTT H— vr — 0 . 


aH m 


H 


(4.21) 


In general, this differential equation cannot be solved analytically, but we can proceed nu¬ 
merically. In any case we need to set up the initial conditions at a high redshift, when dust 
dominates over the dark energy i.e. typically z = 0(100). The initial conditions are, of 
course, model-dependent, but it can be shown that for 2?A-models they reduce to the sit¬ 
uation that comprises the well-known standard ACDM ones for the function and the first 
derivative, viz. 6m{<ii) = o,i and 5'^{ai) = 1, to which we have to add 6'^{ai) = 0. 

Let US recall that the obtained third order differential equation for the matter perturba¬ 
tions does automatically incorporate the leading effect of the DE perturbations and is valid 
under the assumption that the values of the perturbed matter and DE velocity gradients are 
negligibly small. This is the simplest assumption we can make in order to solve the initial 
system of differential equations without introducing additional parameters; and it is certainly 
well justified if we are interested in the physics at scales deep inside the Hubble radius (sub¬ 
horizon scales), i.e k ^ H. From Eq. (4.5) we can see that deviations from this framework 
should scale roughly as ~ We will estimate them in the next section 5.2. 


4.3 Recovering the ACDM perturbations as a particular case 

Let US consider the situation of the ACDM model deep in the matter-dominated epoch when 
we can neglect the A-term. In this case oc a~^ and hence H'{a)/H(a) = —3/(2a) and 
H"{a)/H{a) = 15/(4a^). It follows that Eq. (4.21) boils down in this case to the very simple 
form 

C+l — =<i. ( 4 . 22 ) 

2 a 

whose elementary solution \s 6m = ci a + C 2 + Neglecting the decaying mode and 

setting C 2 = 0 we are led to the standard growing mode solution of the ACDM, i.e. 6 oc a. 
This is a clear indication that the ACDM result is contained in the generalized perturbations 
equation (4.21). 

Alternatively we can prove in a more formal way, in this case using the cosmic time 
variable, that the well-known ACDM linear perturbation equation [46, 47], 

6m + ‘2:H6m ~ ‘^'^Gpm6m = 0 , (4.23) 

or, written only in terms of the Hubble function, 

'6m + “^H^m + H6m = 0 , (4.24) 

is a particular case of (4.20) if we treat the dark energy component as a rigid cosmological 
term with the vacuum constant EoS {ujd = ~1)- First, we perform the time derivative of the 
last equation in order to have a third order one: 

6 m ‘^H6m + 3H6m + H6m = 0. (4.25) 
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To prove that (4.20) and (4.25) are the same when the DE is the traditional A-term, let us 
compute its difference. We find: 


3H6m + 


H5m = 3H 


+ ‘^H6m 



(4.26) 


Upon differentiating the pressure equation (2.2) with respect to the cosmic time under the 
ACDM conditions (viz. = 0) we find H + 3HH = 0, so the difference (4.26) actually reads 

3H (6^ + 2H6^ + H5m) = 0 , (4.27) 

where in the last step (4.24) has been used. This obviously proves our contention. At the 
same time it becomes clear from the proof that if the DE was dynamical (with ud 7 ^ 0) we 
could not have obtained the previous result, which means that the third-order differential 
equation (4.20) is indeed more general than the standard one (4.23) and covers the entire 
class of dynamical DE models with separate (local and covariant) conservation of matter and 
dark energy densities, at fixed gravitational coupling G. 


4.4 Running A and running G with matter conservation 

There is another interesting situation that we would like to mention in passing here in which 
a third order perturbations equation of the same sort is also needed, but under different 
conditions. It was first encountered in Ref. [45] and is characterized by dynamical vacuum 
energy (hence ujd = — 1 ) and conserved matter. This case is, in principle, different from 
the class of models covered in the previous section because the EoS is of the vacuum type. 
However, it differs also in a second aspect, namely in the fact that the gravitational coupling 
G is running or time-evolving along with the vacuum energy density. In this way the Bianchi 
identity can be satisfied since there is a dynamical interplay between the vacuum and G that 
permits the conservation of matter - see [11, 45] for details. We will show now that the 
cosmic perturbations of the matter field for this system are also governed by Eq. (4.20), or 
equivalently by (4.21). In Ref. [45] it was shown that the matter perturbations follow the 
third order equation: 

c + ^(16 - 9D™) + ^i8-an'^ + 3nl - IID^) = 0 , (4.28) 

in which Dm(a) = 8tt G{a) pra{a) / 3H‘^ (a). Equation (4.28) also includes the effects of the per¬ 
turbations in PA) which have been eliminated in favor of the matter perturbations following a 
similar procedure as described in the previous section, except that in this case the dynamics 
of G enters explicitly too. Taking into account that H = aH{a)H'{a) = —47rG(a)pm(a) it fol¬ 
lows that Dm(a) = {—2/3)aH'{a)/H{a)] and using now this relation in (4.28) we immediately 
recover (4.21). This confirms that (4.21) is, in fact, also valid for models with self-conserved 
matter and a dynamical vacuum triggered by the time variation of G. We conjecture that 
(4.21) may also be appropriate for more general models with self-conserved matter and dif¬ 
ferent possibilities for the evolutions of G and the DE, but we will not further pursue the 
framework of G-variable models here - see [21] for a recent analysis. 


- 23 



5 Fitting the PA, VC and PH models to the observational data 

Up to this point we have carried out a thorough study of both the background equations 
as well as of the linear perturbation equations that are obeyed by the class of models under 
consideration. Therefore we are in position to understand the fitting results collected from 
the beginning in Table 1 and that were used throughout to evaluate a number of observables 
associated with the background cosmology in Figs. 1-3. At the same time we need to further 
elaborate on the implications of these dynamical DE models on the structure formation. But 
let us first summarize the methodology employed in our statistical analysis of the cosmological 
data. 


5.1 Global fit observables 

To carry out the fit we have used the expansion history (SNIa+BAO^+BAOrf^) and the 
CMB shift-parameter, as well as the linear growth data, through the following joint likelihood 
function: 

^totip) = C-SNIa X CbAO X CcMB X Cf^g , (5.1) 

in which p = By maximizing the total likelihood or, equivalently, by minimizing 

the joint Xtot function with respect to the elements of p: 

Xtotip) = XsNIa X XbAO ^ XCMB X X/o-g > (5-2) 


we obtain the best-fit values for the various parameters. Concerning the data input, we have 
used the Union 2.1 set of 580 type la supernovae of Suzuki et al. [48]. The corresponding 
XsNia function, to be minimized, is: 


580 

XSNIa(p) = X/ 
i=l 


f^thizi,p) - f^ohsizi) 


(Ti 


(5.3) 


where Zi is the observed redshift for each data point. The fitted quantity /r is the distance 
modulus, defined as/r = m — M = 5 log -|- 25, in which diiz, p) is the luminosity distance: 

dL{z,p) = c{l + z) , (5.4) 

with c the speed of light (now included explicitly in some of these formula for better clarity). 
In equation (5.3), the theoretically calculated distance modulus fith for each point follows 
from using (5.4), in which the Hubble function is the one corresponding to each model, see 
Sect. 3, and we have fixed Hq = 70 km/s/Mpc following the setting used in the Union 2.1 
sample. Finally, ^ohsizi) and (Tj stand for the measured distance modulus and the corre¬ 
sponding Ifj uncertainty for each SNIa data point, respectively. The previous formula (5.4) 
for the luminosity distance applies only for spatially flat universes, which we are assuming 
throughout. 
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We have also made use of the BAO estimators dz(z) and A[z) collected by Blake et al. 
in [50]. The first one can be computed as follows: 


Here 


dz{zi,p) = 

/ OO 

— 


rs{zd) 

Dv{zi) 

cdz 


^(^)t/3 l + lj 


(5.5) 

(5.6) 


is the comoving sound horizon prior to the drag redshift epoch, z^, namely the epoch at which 
baryons are released from the Compton drag of photons and the gravitational stability of the 
former can no longer be avoided by the photon pressure; and 

1/3 


Dviz) = 


(l + z)^Z?l(z)- 


cz 


(5.7) 


Hiz)\ ’ 

is the “dilation scale” introduced by Eisenstein [49]. In this formula, Da{z) = {l+z)~‘^dL{z, p) 
is the angular diameter distance. The other useful BAO estimator is the acoustic parameter, 
A(z), also introduced by Eisenstein as follows [49]: 


A{zi,p) = 




1 


fZi 


dz 


2/3 


(5.8) 


E^l^{zi) [zi Jo E{z)_ 

where Zi is the redshift at which this observable is measured. In all the theoretical expressions 
not related with the SNIa distance modulus we have used the current value of the Hubble 
function given by the Planck Collaboration [3], i.e. EIq = 67.8 km/s/Mpc. The corresponding 
X^-functions for BAO analysis are defined as: 


XBAOdz(P) - ^ 


i=l 


dz,th{Zi,p) dz^ohsi^i 


1 2 


(Jz 


and 


XbaOa 


(p) = E 


6 r 2 

Athiz^i: P) ^obs(^i 


i=l 


^ A,i 


(5.9) 


(5.10) 


where Zi, dz^ohs, crz,i, ^obs and aA,i can be found in Table 3 of [50]. 

The CMB shift-parameter, which is associated with the location of the first peak of 
the CMB temperature perturbation spectrum is given for spatially flat cosmologies by 

dz 





E{z) 


(5.11) 


where z* is the redshift of decoupling. The fitted formulae for the baryon drag redshift, Zd-, 
and the decoupling redshift, z*, are given in [51, 52]. The experimental value for R has been 
taken from [53]. In this way we have constructed the part contributed by the CMB shift 
parameter as follows: 

X^mb(p) - ■ (5-12) 
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Figure 4. Relative differences A(z), Eq. (5.15), between the matter density contrasts <5^(2;) computed 
with the third order differential equation (4.21), which includes DE perturbations, and the result 5m{z) from 
the second order one (5.14), where the DE density enters only at the background level. The curve shows the 
essentially degenerate result valid for all the DA models. 


Finally, we have also taken into acconnt the data on the linear growth rate of clnstering pro¬ 
vided from different sonrces, see specihcally [54-63]. More concretely, we have nsed f(z)as{z) 
and the corresponding x^-function: 

2 / \ f ~ f'^8obs\ 

In recent papers by some of ns [15, 16] we have elaborated a more detailed explanation of all 
these cosmological observables as well as on the htting procednre, and therefore we have left 
these details aside from the present work. We refer the interested reader to the mentioned 
references for more information, see also [10, 11]. In the rest of Sect. 5 we focns on the 
observables nsed to ht the strnctnre formation data, as we feel that they play a very important 
role in onr analysis after we have derived in the previons section the general pertnrbation 
eqnations for the class of self-conserved matter and DE models; and of conrse we report also 
on the impact on onr models. In particnlar, we wish to compare the resnlts of onr ht analysis 
in the sitnation when the DE pertnrbations are inclnded (as analyzed in the previons section) 
and when the DE pertnrbations are not inclnded and the effects of the DE act only on the 
backgronnd history. 

5.2 Structure formation with and without DE perturbations 

As we have explained in Sect. 4.2, we compnte the linear matter pertnrbations dm of the 
models under consideration throngh the third-order eqnation (4.21), which already encodes 
the effect of the DE pertnrbations. When the DE pertnrbations are not inclnded we nse the 
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ACDM form (4.23), which in terms of the scale factor variable can eqnivalently be written as 


0 ^ 


+ 


3 H' 


(5.14) 


with the nnderstanding that the Hnbble fnnction in this expression involves the DE density 
of the corresponding 2?-model in Sect 2. Thns, when nsing this approximation, the DE effects 
enter only at the backgronnd level. It is snrely interesting to compare the resnlts obtained 
in this approximate way with those generated from the more accnrate treatment when the 
DE pertnrbations are compnte from Eq. (4.21). We do not foresee, however, large differences 
between the two treatments since the DE is expected to be mnch smoother than matter, at 
least at the scales where linear strnctnre formation takes place. Still, the fact that we have 
been able to provide a combined treatment of the DE and matter pertnrbations in these 
models gives ns a good opportnnity to test these expectations. 

In Fig. 4 we provide a qnantitative test. We have plotted the relative difference 


A(z) 


- 5m{z) 


5m{z) 


(5.15) 


as a fnnction of the redshift. The density contrast 5^{z) corresponds to the solntion of the 
third order eqnation (4.21) and therefore involves the effect of the DE pertnrbation, whereas 
Smiz) is the solntion to the more conventional second order eqnation (5.14), expressed both 
in terms of the cosmic redshift z = (1 — a)/a. As we can see the differences well inside the 
horizon are very small, of order A ~ 10“®. Let ns recall that the observational data concerning 
the linear regime of the matter power spectrnm lie in the approximate wave nnmber range 
O.Ol/iMpc”^ ^ k < 0.2/iMpc“^, whereas the cnrrent horizon is given by ~ 3000/i“^Mpc 
and hence Hq ~ 3.33 x 10“^/iMpc“^. From onr discnssion in Sect. 4.2 we expect that 
when we take larger and larger scales (i.e. smaller valnes of k as compared to H) potentially 
important contribntions scaling as ~ H‘^/k‘^ can develop. Clearly the ratio sqnared of Hq 
to the minimnm and maximnm valnes of the wave nnmber in the linear regime satisfies 
Hl/k"^ ~ 10-® - 10-3 <C 1 and hence it is natnral to expect that the DE pertnrbations are 
highly snppressed. Onr calcnlation confirms explicitly this fact. Only if we wonld approach 
scales of the order of the horizon snch snppression wonld no longer hold and we conld expect 
sizeable effects from them. 

With this nontrivial test we can say that the DE pertnrbations for the dynamical V- 
models nnder consideration became fnlly nnder control; and in fact we find that at snbhorizon 
scales they do not trigger large deviations from the zeroth order (smooth) DE effects that 
each model already imprints on the backgronnd cosmology. Only at scales comparable to the 
horizon and for non-negligible velocity gradients the DE pertnrbations can play a role. 


5.3 Linear growth and growth index 

In practice it is convenient to investigate the linear strnctnre formation of onr models throngh 
the so-called linear growth rate, namely the logarithmic derivative of the linear density con- 
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trast with respect to the scale factor, /(a) = dln(5m/rfhra [46], or equivalently, in terms of 
the redshift: 


fiz) = -il + z) 


dlndr, 

dz 


(5.16) 


A related quantity is the 7 -index of matter perturbations [46, 65], defined through f{z) ~ 
■ For the ACDM we have 7 ~ 6/11 ~ 0.545, and one typically expects 7 ( 0 ) = 
0.56 ± 0.05 for ACDM-like models [ 66 ]. The growth index can obviously be reexpressed as 


7(^) = 


In f{z) 
In 14 ^( 2 ^) 


(5.17) 


By definition Q,m{z) 
study it reads 


Pm{z) / Pc{z)^ and hence for the class of self-conserved models under 




Pm (1 + zf 

+ zf + Pd{z) 


(5.18) 


in which pd{z) is the DE density of the corresponding 2?-model. Obviously flm(O) = Pm/P^ — 
with p^ = /(SttG) the current critical density. The explicit form for pd{z) 

for the models under consideration has been determined in Sect. 3. 

In recent times, however, a more useful quantity is found to be the weighted growth 
rate of clustering, i.e. f{z)as{z). This quantity has the advantage of being independent of 
the galaxy density bias [57]. Here crs{z) is the rms mass fluctuation amplitude on scales of 
Rs = 8h~^ Mpc at redshift 2 . It can be computed as follows - see e.g. [15] for details: 


asiz) = as, A 


dmiz) 

dm 


e+nsT^n^,k)W^(kRs) 

k'^+^^T^{n^,k)W^{kRs) 


-,1/2 


(5.19) 


with Us = 0.968 ± 0.006 and = 0.815 ± 0.009 (cf. [3]). The smoothing function 
W{kR) = 3{smkR — kRcoskR)/(kR)^ is the Fourier transform of the following geomet¬ 
ric top-hat function with spherical symmetry: /top hat (^') = 3/{ATrR^) 9 h{ 1 — r/R), where 
9h is the Heaviside function. It contains on average a mass M within a comoving radius 
R = {3M/Att pmY^^■ In the above expression T(Hm\ k) is the BBKS transfer function [47, 67]. 
Introducing the dimensionless variable x = k/k^q^ in which keq = aeqH{aeq) is the value of the 
wavenumber at the equality scale of matter and radiation, we can write the transfer function 
as follows: 


T{x) = 


ln(l -h0.171x) 
0.171a; 


1 -F 0.284X -F (1.18x)2 -k (0.399x)^ -F (0.490x)^ 


-1/4 


(5.20) 


Notice that keq is a model-dependent quantity that must be derived for each of the models 
we are studying in this paper. For the ACDM we have the standard result [47] 


(ACDM) : 


ki\ = Hnn)LJ-^ 


•'eq 



2 -no^b 


(5.21) 
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Figure 5. Left: Comparison of the observed and theoretical evolution of the weighted growth rate/( 2 ;)ct 8 («) 
for the ACDM versus the VA and DC2 models; Right: The evolution of the growth index (5.17) for the same 
models in a wide span of 2 so as to display the vertical asymptotes at the points where the DE density pn ( 2 ) 
vanishes in each case and the 7 -index (5.17) becomes singular. The bar in the models correspond to using the 
barred fitting parameters of Table 1, whose meaning is explained there. 


However, for the DA and DCl we have to introduce a correction factor, which is analytically 
calculable. We hnd: 


{DA): 


(PCI) : 


keg = k, 


3/2 


+ 


1/2 


-|l/2 


4a + 3{l — iy) 1 + a — u 


k — 

^eq — 


eg 




(5.22) 

(5.23) 


As we can see, the correction factor in (5.22) depends separately from v and a, not just on the 
difference. The corresponding formulas for DAl and DAS models is found by setting a = 0 
and 1 ^ = 0, respectively, in the hrst expression. Notice that the hrst formula above reduces 
to the second for a = 0, as expected from the fact that at high z the model PCI behaves as 
PAl. This also explains why the effect of the e-term for PCI is negligible in this regime and 
does not appear in (5.23). As for the correction factor for model PC2 it is derived also from 
(5.22). Recall that once is hxed from the ht value, the parameter a becomes determined 
for this model under the conditions discussed at the end of Sect. 3.2. 

The model-dependent corrections to keg do modify the shape of the transfer function 
and be of importance in some cases. It can be instructive for the reader to compare the above 
correction formulas with the ones determined for the corresponding dynamical vacuum models 
in [16] - see Eq. (59) in that paper. One may be a bit surprised at hrst that the formula for 
the correction factor found there for Cl (denoted in that reference as model II) is considerably 
more complicated than the one found here for PCI. The reason is that the matter density for 
the former model has, at variance with the latter, an anomalous behavior at high redshift - cf. 
our discussion in Sect. 3.2, particularly Eq. (3.28). Even so the correction indicated in (5.23) 
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for PCI can be significant because for this model v is not a priori a negligible parameter (cf. 
Table 1). As a matter of fact, even for models of the PA-subclass the correction to k^q is of 
some signihcance for the overall ht, e.g. it diminishes slightly the value with respect to 
the situation without correction. 

5.4 Results 

Let us now further discuss the specihc results obtained for the P-models under study. Table 
1 , repeatedly referred to in the previous sections, collects in a nutshell the main numerical 
output describing the outcome of our analysis. The basic traits of the background history of 
these models have been presented in Sect. 3. Here we mainly focus on the structure formation 
results. They are represented in Figs. 5-7. Furthermore, the contour plots combining all the 
observational data (including of course the structure formation data extracted from the linear 
growth rate and associated quantities) is indicated in Fig. 8 . Next we discuss these hgures 
in turn. 

In Fig. 5 (left) we consider models PA and PC2, together with the concordance model 
ACDM. We plot their theoretical prediction of f{z)as{z) versus the data points (provided 
in the above mentioned references [54-63]) as a function of the redshift and for the best ht 
values of the parameters in Table 1. We use the barred parameters of the table, as in this 
case we take into account the structure formation data in the ht. Notice that since model 
PA2 interpolates between PAl and PA3 we have plotted the last two only in order to show 
the range that is covered. The htting results for PA2 recorded in Table 1 correspond to 
hxing a = —ly > 0. In fact, since a must be positive (cf. Sect. 3.1) the previous setting is 
illustrative of how to break degeneracies among the parameters in a way that is consistent 
with the sign of a. Notice that, at this point, this is necessary since the v and a dependence 
in e.g. equations (5.22) and (5.23) is no longer of the form v — a. Other choices compatible 
with a > 0 give very similar results. 

The reason why we have included PC2 also in Fig. 5 stems from our discussion in 
Sect. 3.2, where we showed that this model can mimic the ACDM near our time, so it is 
natural to joint it in the same group of models approaching the concordance cosmology now. 
We can see that the three curves thrive quite well and stay together slightly below the ACDM 
line. They actually provide a better ht than the ACDM (cf. Table 1). We will further qualify 
the magnitude of this improvement in the next section. On the other hand, in Fig. 5 (right) 
we depict the growth index, 7 ( 2 ;), dehned above. Here the three models are seen to approach 
together only when z —?■ 0, i.e. around the current time. While PAl and PA3 still remain 
quasi-glued at all points, PC2 departs signihcantly from them at larger redshifts. So the 
mimicking of the ACDM by PC2 is actually not better than that of the PA-models. Let us 
also note the presence of vertical asymptotes, which are connected with the vanishing of the 
denominator in Eq. (5.17). These are the points where the DE density pd{z) becomes zero 
(confer Fig. 1), and hence t 1 from Eq. (5.18). Recall also from Fig. 2 that at these 
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Figure 6. Left: The evolution of the linear growth rate of clustering, Eq. (5.16), specifically for the 
general DC 1-model and the DH one (i.e. the linear model po ~ H) both for the barred and unbarred fitting 
parameters indicated in Table 1. For reference, the ACDM model is also included; Right: The corresponding 
curves for the weighted growth rate f{z)as{z). The dash-dotted line (in blue) in both plots corresponds to 
the best-fit values of the parameters for the situation when the structure formation data are not used. This 
plot makes apparent the delicate situation of the DCl-model when the unbarred fit parameters are used, 
particularly in terms of f{z)crs{z) rather than just in terms of f{z) - see also the text. 


same points the EoS parameter of the DE becomes singular and develops a corresponding 
asymptote. 

In Fig. 6 we compare the situation of the DCl and DH models (together with the ACDM 
as a reference). We have separated these P-models from the rest because we know that they 
behave differently. Let us first focus on the plot on the left, where we display the linear growth 
function f(z) for them, Eq. (5.16). In the case of DCl we include both the curve using the 
barred parameters and the unbarred ones (cf. Table 1). We can see that in the latter case the 
"DCl curve is displaced only slightly below the others. However, when we inspect the figure 
on the right, corresponding to the weighted function f{z)as{z), i.e. the growth function 
appropriately rescaled with the rms mass fluctuation amplitude at each redshift, the same 
DCl curve now strays downwards quite significantly from the rest. As it turns, the weighted 
linear growth f{z)as{z) appears specially sensitive to that. It demonstrates that when we 
attempt to fit DCl without using the structure formation data the model immediately gets 
in tension with them, quite in contrast with the situation with the DA models. 

The anomalous behavior of DCl (and, a fortiori, of DYl) is particularly clear from 
the numbers in Table 1, where the difference between the unbarred and barred values of 
the fitting parameter (viz. u ~ —0.64 and i> ~ —0.35) is abnormally large, in contrast to 
the other models. Such behavior is also reflected in the large increase of recorded in 
Table 1, for this model when we compare x^, that is to say, the “reduced x^ value” (the 
value which is computed without including the linear growth X/o-g contribution) with the 
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Figure 7. Magnified view of the growth index plot (cf. Fig. 5, right) for the local range 0 < z < 2. We have 
superimposed also the growth index for the general I?Cl-model (dash-dotted blue line) and the linear model 
VH, dash-dotted grey line). The curves have been obtained under the same conditions of the two previous 
figures. 


unreduced one ~ when the growth data are counted in the computation. The latter can be 
computed either without optimizing the ht to the growth data (yielding /dof = 880.74/600) 
or after optimizing it (rendering x^/dof = 635.23/600). In both cases the result is much 
larger than Xr/dof = 563.86/584. Obviously the second result is better than the hrst, but 
both are extremely poor since they have to be compared with the respective ACDM values, 
X^/dof = 584.91/608 and x^/dof = 584.38/608. This is in contrast to the situation with the 
VA models, where we can check that the difference between the reduced Xr value and the 
corresponding x^ values computed with or without optimizing the fit to the growth data is 
by no means as pronounced as in the PCI and PH cases. 

We can actually retrace the same features described above if we look at Fig. 7. This is a 
magnified view of the growth index plot in Fig. 5 (right) for the local range 0 < z < 2, where 
we have superimposed also the growth index for the general PCI-model and the linear model 
PH. These lines are in correspondence with the ones in Fig. 6 and highlight once more the 
anomalous behavior of the PCI-model and its exceeding sensitivity to the structure formation 
data. The current growth index value 7 ( 0 ) > 0.66 (resp. 7 ( 0 ) > 0.72) that one can read off 
Fig. 7 for PCI when the structure formation data are used (resp. not used) in the fit, is 
clearly too large to be acceptable. The vertical asymptote at z ~ 1.5 (where Pd{z) vanishes) 
is much closer here because the best-fit value for the unbarred v is much bigger in absolute 
value than i>. As for PH it has no asymptote because the corresponding Pd{z) never vanishes 
(cf. Fig. 1). This model also separates from below the ACDM line. Models PA, instead, 
perform quite well and meet sufficiently close the ACDM horizontal line at 7 ~ 6/11 ~ 0.545 
for z —)■ 0 . 


- 32 - 









In short, the unfavorable verdict for DCl seems to have no cure: when we enforce 
the adjustment of the growth points, then it is the expansion history data that becomes 
maladjusted; and conversely, if we strive to adjust the expansion history data, then it is the 
structure formation data that get strayed. In both cases the value of "DCl becomes too 
large in comparison to the other models. The exceeding sensitivity to the growth rate and 
associated quantities is symptomatic of its delicate status as a candidate model to describe 
the overall observations. And, of course, everything we have said for "DCl applies to PH, 
which is the particular case = 0 of the latter. The model PH has no free parameter in 
the vacuum sector and therefore there is nothing we can do to improve its delicate situation, 
which in the light of Table 1 appears to be the less favorable one among the models under 
study. 

From the plots in Figs 5-7 and the statistical results of Table 1 we can say that the 
PA models are capable of describing the overall set of data in a way which is perfectly 
competitive with the ACDM. Model PAl, in particular, represents the simplest realization 
with one parameter, u. The same can be said of model PA3 with the parameter a. These 
models are similar because the time evolutions induced by and H are comparable. As 
for model PA2 it is the most general of the PA-subclass and interpolates between the two 
previous cases. All of them offer a better fit quality than the ACDM. In the next section we 
further qualify this assertion. 

5.5 Discussion 

Let us now focus on Fig. 8, where we display a detailed representation of the fitting procedure 
we have followed in order to pin down the most favorable region of the various DE models 
under study. In that figure we specifically consider the case of model PAl, where we can 
carefully appraise the way we have combined all the data on expansion history and structure 
formation mentioned in Table 1. Thanks to the juxtaposition of the permitted regions defined 
by each observable we find that a well delimited domain emerges. This bounded, elliptically 
shaped, area is what we may call the physical region for each model. The upper plot on the 
left of that figure displays all the lines involved in the projection of the final physical region; 
right next we show a plot with the boundary lines corresponding to the CMB shift parameter 
and the weighted growth rate f(z)as{z)‘, then the lower plot on the left highlights the edges 
associated with the two types of BAO data involved, together with the borders of the area 
allowed by the supernovae (SNIa) observations; finally, the net intersection of all of them is 
represented by the lower plot on the right, in which we depict the la, 2a and 3a domains 
within the physical region. The corresponding plots for models PA2 and PA3 are very similar 
and are not shown. 

A most remarkable feat emerging from our analysis is the fact that the PA-subclass of 
dynamical DE models is capable of improving the fit quality of the ACDM. This observation 
is specially significant if we bare in mind that the physical region in Fig. 8 is mostly located 
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Figure 8. Likelihood contours in the plane (for the values —2\nC/Crnax = 2.30, 6.16,11.81, 

corresponding to Icr, 2<j and 3cr confidence levels) for the HAl-model using the full expansion history and 
CMB shift parameter data (BAO+SNIa+CMB). Models VA2 and VA3 present very similar plots, as can be 
expected from the statistical analysis of Table 1, and for this reason we have denoted VeS the vertical axis of 
the plot. In all cases the iZeft = 0 region (ACDM) is disfavored at ~ 3cr level. 


outside the domain of the concordance model, namely it lies roughly ~ 3cr away from the 
z^eff = 0 line representing the ACDM; specifically, we can estimate with accuracy that 84.27% 
of the area of the 3cr domain is below the = 0 line. Models DA, therefore, represent a 
challenge to the ACDM since they are phenomenologically quite competitive. Somehow they 
“break the ice” inherent to the rigid character of the cosmological A-term, showing that the 
possibility of a dynamical vacuum is not only more natural from the theoretical point of view 
but also more instrumental for adjusting the cosmological observations. 

In order to assess the statistical quality of our fits, and the competition score between the 
models, we have displayed their global value per degree of freedom in Table 1. However, 
in the same table we provide also (in the last columns) the value of another related statistics, 
which we now describe. It is the so-called Akaike Information Criterion (AIC) [36], which is 
particularly useful for comparing competing models describing the same data and is amply 
used in many areas of science — see e.g. the book [68]. We wish to use the AIC here to 
compare the D-models with the ACDM. Such criterion is formulated directly in terms of the 
maximum of the likelihood function, £, and is defined as follows. Let Up be the number of 
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independent estimable parameters in a given model and N the sample size of data points 
entering the fit. If £max is the maximnm valne of the likelihood fnnction, the corresponding 
AIC valne is defined as follows: 

AIC = -2In+ 2„, + , (5.24) 

The comparison criterion based on this statistics is the following: given two competing models 
describing the same data, the model that does better is the one with smaller AIC valne. Notice 
that there is a kind of tradeoff between the first term on the r.h.s. of (5.24) and the other 
two. The first term (the one carrying the maximnm likelihood valne) tends to decrease as 
more parameters are added to the approximating model (since Tmax becomes larger), while 
the second and third terms increase with the nnmber of parameters and therefore represent 
the penalties applied to onr modeling when we add more and more parameters. In particnlar, 
the second term (2np) represents an universal penalty related only to the increase in the total 
nnmber of independent parameters, whereas the third term is an extra penalty applied when 
the sample is not snfhciently large as compared to the nnmber of independent parameters. 
For large samples N ^ Up (typically for N/up > 40) the third term becomes negligible and 
the above formnla simplifies to 


AIC =-21n£max + 2np. (5.25) 

This is actnally the sitnation that holds good in onr case. Indeed, for the ACDM we have one 
independent parameter (up = 1 ) represented by whereas for PAl, for instance, we have 
Up = 2, namely (fl^, v). In the case of 'DA2 we have (12^) hence np = 3, bnt since 

we fix a relation between v and a. we have again np = 2. On the other hand we can see from 
Table 1 that N ~ 590 in all cases, and therefore the condition A^/rip > 40 is amply satisfied 
and the nse of the more simplified formnla (5.25) is fnlly jnstified in this case. For Ganssian 
errors, the maximnm of the likelihood fnnction can be expressed in terms of the minimnm of 
and therefore Eq. (5.25) can be reexpressed as 

AIC = xLn + ^np . (5.26) 

In practice, to test the effectiveness of models Mi and Mj, one considers the pairwise difference 
(AIC increment) (AAIC)ij = (AIC), — (AIC)j. The larger the valne of Aij = |A(AIC),j|, the 
higher the evidence against the model with larger valne of AIC, with Ajj > 2 indicating a 
positive snch evidence and Ajj > 6 denoting significant snch evidence. Finally the evidence 
ratio (ER) against one model whose AIC valne is larger than another is jndged by the relative 
likelihood of model pairs, Tj/Tj, or eqnivalently by the ratio of Akaike weights Wi = e(AiC)i/ 2 ^ 
Thns, given two models Mi and Mj whose AIC increment is Ajj, the evidence ratio (ER) 
against the model whose AIC valne is larger is compnted from ER = Wijwj = . 

From Table 1 we can derive the AIC increments (AAIC),j and evidence ratios when we 
compare the fit qnality of models i ='DA. ,PC ,2?!! with that of j = A (ACDM). We find that 
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for all PA models (AAIC)jA ^ 9 against the ACDM. It follows that when we compare any 
PA model with the ACDM the typical evidence ratio against the latter is typically of order 
ER~ 90. Worth noticing is also the result of the fit when we exclude the growth data from 
the fitting procedure but still add their contribution to the total ~ h corresponds to the 
unbarred parameters in Table 1. This fit is of course less optimized, but allows us to risk 
a prediction for the linear growth and hence to test the level of agreement with these data 
points. It is noteworthy that in the case of the PA models the corresponding AIC pairwise 
differences with the ACDM remain similar to the outcome when the growth data enters the 
fit. Therefore, the ACDM appears significantly disfavored versus the PA-models in both cases 
(i.e. with barred and unbarred parameters in Table 1). 

Quite another story is the situation with the PCI-model, which does not seem capable to 
adjust simultaneously the expansion history observables and the structure formation data. A 
glance at Table 1 shows that when we compare it with the ACDM, with barred parameters, we 
find (AAIC)ia ^ 53 against PCI, whereas with unbarred ones we get the even worse result 
(AAIC)ia ^ 296 . The respective evidence ratios against PCI versus the ACDM under these 
two circumstances are of course extremely large, the smallest one being ER rsj 10^^. Thus, 
quite obviously, from the point of view of the Akaike Information Criterion the situation of 
the PCI model is desperate, and that of PH is even more dramatic. Not so for the PC2 
model, but we will assess its viability more carefully at the end. 

The upshot after making use of the AIC is that the competitive position of the PA 
models became even more prominent from the point of view of phenomenology, whilst the 
status of PCI and PH against observations appear essentially as terminal. This fact may 
have nontrivial implications for cosmological model building. For instance, the PCI model 
has been studied in some formulations relating QCD with cosmology and it goes sometimes 
under the name of “QCD ghost dark energy model” [33-35]. Although some formulations of 
these models are certainly of theoretical interest and could possibly be modified to become 
in better harmony with observations, the strict phenomenological analysis of models of DE 
containing a linear power of the Hubble function and a quadratic power of it, with no additive 
constant term {Cq = 0) - i.e. the strict PCl-form indicated in our model list (2.10) — leads 
to the inescapable conclusion that they are unable to account for the observations. The 
particular case in which the DE contains only the linear term in H, i.e. po ~ H, suffers of 
an even more hapless fate. Such linear form was first proposed in Ref. [32] and subsequently 
was adopted and adapted to different purposes by different authors, until it eventually gave 
rise to the extended PCI form. Incidentally, the two PC models (2.10) and the linear (2.11) 
were previously shown to be phenomenologically problematic as well when treated as vacuum 
models, see [15] and in particular [16]. Here we have shown that the phenomenological 
conflict also persists when they are treated as P-models with self-conserved dynamical DE. 
If we compare Fig. 5 of Ref. [15] with the present Fig. 6 we can see that the PCI-model 
actually does better than their vacuum counterpart, the Cl-model (whose lack of power for 
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structure formation near our time is quite dramatic, as shown in that reference). Still, the 
overall performance of the DC 1-model remains helpless in front of the ACDM and its DA 
companions. 

We would like to point out here that our assessment about the DCl and DH models 
is definitely much more pessimistic as compared to the analysis presented e.g. in [35] - see 
also the review [64]. Take e.g. DCl; the basic difference that we have been able to trace 
with respect to these authors is that we do not find that such model can pass the structure 
formation test when DCl is adjusted to the background data, and vice versa. We suspect 
that in the aforementioned references the linear growth data points have directly been fitted 
with the DCl model without perhaps verifying if the density perturbations of the model were 
able to meet the ACDM behavior at large redshift, namely 5 (x a. In our case we fitted 
the data after imposing such boundary condition, and we found that the model failed (by 
far) to reproduce the observations, as we have reported in detail throughout our work. The 
results exhibited in Figures 6-7, combined with the comparatively poor statistical output of 
the DC 1-model recorded in Table 1, altogether make a quite eloquent case against the delicate 
health of that model in front of observations. The situation with DH, which is a particular 
case of DCl, is even worse since it has one less degree of freedom to maneuver. In short, from 
the present comprehensive study we conclude that the DCl and DH models are ruled out. 
This was already the firmly inferred conclusion when these models were treated within the 
vacuum class [14-16],^ and here we can confirm that it is also the same unfavorable situation 
in the context of the D-class. 

Let us finally comment on the viability of the DC2-subclass, which was shown to fit the 
data reasonably well (cf. Table 1). It is sometimes associated with the so-called entropic-force 
models [31]. We have seen in Sect. 3.2 that it can provide a reasonable fit to the low redshift 
data (including linear growth) and therefore performs better than the DCl and DH models. 
However, this conclusion cannot be placed out of context, meaning that the structure of the 
DC2 model - as in fact of any other one in the list of DE models given in (2.9-2.11) - is 
not sufficient to conclude whether the model is viable or excluded; for it also depends on 
whether there is exchange of energy with matter or not. Thus, despite the DC2 model (as a 
self-conserved DE model) provides a reasonable fit to the low-z data in Table 1, its vacuum 
counterpart - i.e. the model C2 with ojd = —1 and matter non-conservation studied in detail 
in [14, 15] - was shown to be ruled out, already at the background level, even though it has 
exactly the same structure as a function of H. 

Unfortunately DC2 has ultimately, too, a very serious drawback as self-conserved DE 
model, which we have preliminary announced at the end of Sect. 3.2. While DC2 tends to 
mimic the ACDM for a while, which is of course a positive attribute, such feature is limited 
to a period around the current epoch and cannot be extended to the radiation-dominated 

^See in particular the thorough analysis of Ref. [16], including the additional details provided in the 
Appendix of that paper and references therein. 
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era. The reason can be seen on inspecting the radiation term in Eq. (3.7), namely the one 
proportional to ^1^/(1 — v + 4q:/3). This term does apply to PC2 as well (with Cq = 0), 
and the problem appears becanse the \ow-z data natnrally choose ~ a ~ 1 (cf. Table 1) 
so as to mimic as mnch as possible the ACDM. However, this is only possible at the price 
of “renormalizing” the size of the radiation coefficient from to 3H(!/4. Obvionsly this 
is a rather significant modihcation, in which 25% of the normal content of the radiation is 
missing. Thns, deep in the radiation-dominated epoch, the model becomes anomalons and one 
can foresee a signihcant departure from the ACDM. For this reason and despite its ostensible 
snccess at low energies we do not place this model in the list of onr favorite P-models. A 
more detailed analysis (which we do not deem necessary here) wonld reqnire to appropriately 
modify the standard formnlas existing in the literatnre e.g. for compnting the deconpling and 
baryon drag redshifts, as well as the corresponding modihcation of the BAO^^ observable (cf. 
the discnssion of these formnlas in [15] and references therein). 

We close this section by mentioning the analysis of Ref. [69] on a similar model as onr 
2502 . The anthors do not seem to have detected the above mentioned important problem, as 
in fact they did not compnte the radiation contribntion. In addition, these anthors incorrectly 
compare their resnlts with those of [14] withont realizing, or at least clarifying, that the two 
models are very different; namely, in one case it is a Cl-type vacnnm model [15] in interaction 
with matter, whilst in the other it is a PCI-model with self-conserved matter and DE. Be 
as it may, in the light of the resnlts presented here it becomes clear that all the models with 
Cq = 0, whether in the vacnnm class or in the P-class, are actnally exclnded and cannot be 
nsed to snpport an alternative formnlation of holographic DE . 

6 Conclusions 

In this work we have thoronghly analyzed the P-class of self-conserved dynamical dark energy 
models. These models are based on a dynamical DE density pD which takes the form of a 
series of powers of the Hnbble rate, H, and its derivatives. To be consistent with the general 
covariance of the effective action of qnantnm held theory in cnrved spacetime, we can expect 
that the series of powers mnst involve an even nnmber of cosmic time derivatives of the scale 
factor only. Thns, if we shonld apply strictly this recipe for the cnrrent Universe, only the 
hrst powers np to H and wonld be involved in the strnctnre oi pD, inclnding or not a 
constant additive term Cq. However, in order to cover a wider variety of possibilities that 
have also been addressed in the literatnre from different theoretical perspectives, we have 
admitted also a snbclass of models in which the linear term in H is inclnded for Cq = 0. This 
greater hexibility in the general form oi pD has allowed ns to compare the performance of this 
restricted set of models with respect to those that are theoretically most preferred, namely 
the Co / 0 ones. 

An additional characteristic of the P-class is that pD is locally and covariantly self- 
conserved. At hxed valne of the gravitational constant, this entails the simnltaneons local 
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covariant conservation of matter as a necessary condition to satisfy the Bianchi identity. The 
"D-class models are thus enforced to have a dynamical equation of state with a nontrivial 
evolution with the cosmic expansion, cod = ujd{H), and therefore it is different from the 
previously studied dynamical vacuum class [15] in which = —1 and the vacuum exchanges 
energy with matter. The two dynamical DE types share the same formal functional depen¬ 
dence of pd on H, but the behaviors are substantially different, already at the background 
level, as attested by the analytical and numerical solutions in each case. 

We have also studied for the first time (to the best of our knowledge) the generic set of 
linear matter and DE perturbations for a system of self-conserved cosmic components, and 
studied the conditions by which the system can be transformed into an equivalent third order 
differential equation for the matter perturbations. We have demonstrated that the situation 
of a rigid A-term, i.e. the ACDM concordance model, appears as a particular case of that 
general framework. We have subsequently applied it to the entire P-class of DE models and 
upon solving the general perturbation equations we have compared the solution with the 
corresponding results obtained when the DE enters only at the background level. In this way 
we have been able to have good control on the reach of the DE perturbations and explicitly 
confirmed that they play a small role when the considered scales are well under the horizon. 

After a detailed study of the background history and the cosmic perturbations of the var¬ 
ious models in the P-class, we have been able to identify the subclass of the DA-models as a 
most promising one insofar as it provides an excellent fit to the overall data on Hubble expan¬ 
sion and structure formation. These are the theoretically favored P-models for which Cq 7^ 0 
- cf. Eq. (2.9). The fit quality rendered by them has been shown to be significantly better 
than that of the ACDM. Most conspicuously, using the Akaike Information Criterion [36] as 
a method to compare competing models describing the same data, we find that the evidence 
ratio in favor of the PA-subclass (namely the relative likelihood of these models as compared 
to the concordance model) is of order of a hundred. We also find that the physical region of 
the parameter space for the PA-models lies ~ 3a away from the ACDM region, and therefore 
the two models can be clearly distinguished. 

Using the same testing tools we have reached the firm conclusion that all of the P-models 
with Co = 0 are strongly disfavored, in particular the linear model pD ~ H. Furthermore, 
among the theoretical models existing in the literature that become automatically excluded 
by our analysis we have the so-called entropic-force models and the QCD-ghost formulations 
of the DE (cf. Sect. 5.5 for details and references). 

At the end of the day the most distinguished dynamical 2?-models, both theoretically 
and phenomenologically, are those in the 2?A-subclass ~ viz. the set of models which endow 
the DE of a mild dynamical character and at the same time have a well-defined ACDM limit. 
These models improve significantly the fit quality of the ACDM, showing that a moderate 
dynamical DE behavior is better than having a rigid A-term for the entire cosmic history. 

We have found that the favored PA-subclass has also the ability to mimic the quintessence 
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behavior and it could even provide a possible explanation for the phantom character of the 
DE at present, as suggested by the persistent region projected below the w/j = — 1 line (the 
so-called phantom divide) from the fits to the recent and past cosmological data. There is of 
course no significant evidence of this phantom character, as the physical region includes the 
ujd = —1 line. However, if in the future more accurate observations would insist on singling 
out the domain below the phantom divide, then the dynamical DE models in our "DA-subclass 
could provide a simple explanation without need of invoking true phantom fields, which are 
of course abhorred in QFT. 

Let us conclude by emphasizing our main message. The dynamical DE models treat 
the DE density as a cosmic variable on equal footing to the matter density. In a context of 
an expanding universe this option may be seen as more reasonable than just postulating an 
everlasting and rigid cosmological term for the full cosmic history. Furthermore, the subclass 
of "DA-models favored by our analysis furnishes a theoretically consistent and phenomenolog¬ 
ically competitive perspective for describing the dark energy of our Universe as a dynamical 
quantity evolving with the cosmic expansion. The structure of the DE density in these mod¬ 
els as a series of powers of the Hubble rate and its time derivatives is suggestive of a close 
connection with the QFT formulation in curved spacetime. We hope that they will help to 
better understand the origin of the cosmological term from a more fundamental perspective, 
and hopefully they might eventually shed some light as to the nature of the cosmological 
vacuum energy and its relation with the quantum vacuum of modern gauge field theories. 
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